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Abstract 

We study, in the framework of the Boltzmann-Nordheim equation (BNE), 
the kinetic properties of a boson gas above the Bose-Einstein transition tem- 
perature T c . The BNE is solved numerically within a new algorithm, that has 
been tested with exact analytical results for the collision rate of an homoge- 
neous system in thermal equilibrium. In the classical regime (T > 6 T c ), the 
relaxation time of a quadrupolar deformation in momentum space is propor- 
tional to the mean free collision time T re i ax ~ T -1 / 2 . Approaching the critical 
temperature (T c < T < 2.7 T c ), quantum statistic effects in BNE become 
dominant, and the collision rate increases dramatically. Nevertheless, this 
does not affect the relaxation properties of the gas that depend only on the 
spontaneous collision term in BNE. The relaxation time r re i ax is proportional 
to (T — T c )~ 1 / 2 , exhibiting a critical slowing down. These phenomena can be 
experimentally confirmed looking at the damping properties of collective mo- 
tions induced on trapped atoms. The possibility to observe a transition from 

collisionless (zero-sound) to hydrodynamic (first-sound) is finally discussed. 
PACS numbers: 05.30.Jp, 05.20.Dd, 51.10.+y 



1 



Typeset using REVTpjX 



The study of the kinetic properties of a (non-condensate) Bose-Einstein (B-E) gas, being 
interesting in itself, is becoming particularly important in relation with the recent experi- 
mental observations of the Bose-Einstein condensation (BEC) of trapped atoms and 
excitons ||. In the experimental setups of a gas of alkali atoms is trapped and cooled 
down below the BEC critical temperature (T ~ 100 nK for iV ~ 10 4 atoms in jlj]). The 
cooling is achieved tuning the trapping magnetic field in order to allow the most ener- 
getic particles to escape, with relaxation processes(essentially one- and two-body collisions) 
driving the system to new equilibrium states at lower temperatures. It is clear that the 
understanding of the dynamics of the condensation, requires a kinetic model describing the 
interplay between the time scale associated with the evaporative cooling, the relaxation pro- 
cesses and the emerging of the condensate f5H8|. Such a complete analysis is still lacking 



One possible scenario is to assume that the time scales associated with the growth of the 
condensate are large compared with the relaxation to equilibrium. The condensation will 
have a quasi-static character allowing the description of the system in the framework of the 
grand-canonical statistical ensemble P-|TT|]. On the other hand, for such dilute systems (the 
mean inter-atomic distance is much greater then the scattering length), the equilibration 
times can be far longer than the time associated with the growth of the condensate. In this 
case, non-equilibrium processes would be crucial in understanding the BEC. Yet, the finite 
life-time of the trapped atoms is in competion with the time necessary for the emergence of 
the condensate. This particular problem explains the so far unfruitful search of condensation 
in a trapped gas of Hydrogen atoms fll^fl . A further intriguing possibility, is the observation 
of the transition from collisionless (zero-sound) to hydrodynamic (first-sound) modes, that 
occurs when the number of collisions is enough to reach local thermodynamic equilibrium. 

Several theoretical papers address the dynamics of a weakly interacting B-E gas. Levich 
and Yakhot consider its evolution in the presence of a thermal bath ||13|| , finding that the 
time for condensation is infinite. Snoke and Wolfe [TJ[] solve numerically the BNE for an 



homogeneous, thermally isolated, boson gas. They demonstrate that in the classical regime 



(high temperatures and/or low densities) the Maxwell-Boltzmann distribution is achieved 
within few inter-atom scattering times. At sufficiently high densities, in the degenerate 
regime, the collision rate increases, but the number of collisions necessary to relax the 
system increases as well. The net result is that the relaxation time remains proportional to 
the classical mean collision time. Again, the appearance of a B-E condensate (an infinite 
value of the distribution at zero momentum) is not achieved. However, a different conclusion 
is reached by Semikoz and Tkachev []15| |, who found the signature of the B-E condensation 



in a finite time. 

The authors of consider an isotropic deformation in momentum space. In this case 

it is possible to simplify greatly the collision integral in BNE, making feasible the numerical 
calculations in a grid of discrete energy points. 

In this paper we study, as a function of the temperature, both the collision rate of a 
B-E gas at equilibrium and the relaxation time of a quadrupolar deformation in the mo- 



mentum space. Differently from the cases studied in [T^|T^] we need to solve the BNE 
in the full momentum space, without assuming the isotropy of the system. To this end, 
we develop a new numerical algorithm, that we test with analytical results. Our analsys 
indicates that for a quadrupolar (non-isotropic) deformation the relaxation times diverges 
as r re i ax oc (T — T c ) -1 / 2 , exhibiting a critical slowing down. This could seem in contradic- 
tion with the dramatic increase of the collision rate that occurs approaching the critical 
temperature. However, the effect is easily understood as due to the cancellations between 
the stimulated "gain" and "loss" collision terms in the BNE, with the relaxation time that 
depends only on the "classical" spontaneous collision term. As a consequence, close to the 
critical temperature T < 2.7 T c the system exhibits zero-sound modes, and the simple crite- 
ria given in the literature that state the hidrodynamic limit occurs when ujt co u « 1, with 
u> the frequency of the collective motions, clearly breaks down. 

We note that the BNE approach cannot describe the onset of the B-E phase transition. 
In fact, the emerging of coherence (with all atoms in the condensate having a common, well 
defined phase) requires higher order correlations @ not included in BNE, that is derived in 



the random phase (semiclassical) approximation ||. 

The experimental study of kinetic properties of the B-E gas is quite active. Collective 
motions of a trapped condensate have been induced using time- dependent magnetic fields 
with different multipolarities []I6 , |T7| . The frequencies fl8|-|20|| and the relaxation times [B 



of these resonances have been studied in detail, and the interest is now focusing on their 



finite-temperature properties [^T 



The modification of the classical Boltzmann equation in order to take into account the 



effect of Bose statistics in the collision integral was first suggested by Nordheim [p2 
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with g the spin degeneracy that, later on, will be set g = 1. The collision rate W is related 
with the scattering cross-section as 

W=**£ (2) 

where a = 8ira 2 , with a the scattering length. Here we develop a Boson Monte Carlo 
Dynamics (BMD) method to solve this equation for the homogeneous case (that can be 
generalized to solve the BNE eq.(l) in the full phase-space). We parametrize the one-body 
distribution function as a sum of "weighted-test-particles" (WTP): 

N TP 

n(p,t) = c W (p-Pi( t )) 

(2nhf (3) 
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where p is the density of the system, and N TP is the total number of "test-particles". We 
choose the weight functions w(p — pi) = 5{p — pi) where pi{i) is the time dependent position 
in momentum space of the test-particle i. 

The time evolution of the particles swarm proceeds along the following steps: 



1) the test-particles are randomly distributed in momentum space according to the initial 
distribution n(p, t — 0) 

2) at each time-step At, we randomly choose two test-particles and we define a 
collision probability as: 

P coll = 1 - exp[- ^ Pi ~ Pi ^ p a FtfrfJ) At] (4) 

where F — [1 + n(p-)][l + 

3) A random number R = {0, 1} is compared with P co u. if R < P co « the collision is 
accepted. 

4) The new momenta and p- that the two particles would have if the collision were 
accepted are randomly chosen imposing the conservation of energy and momentum. 

4) Repeat the process from 2) until all test particles have had the chance to collide in 
that time step. 

The total energy, momentum and number of particles are exactly conserved at all time- 
steps. It is clear that At (~ 0.01 s) and N tp (~ 1000) should be chosen to reach the 
convergence of the results. It is especially important for the time step At to be small 
enough in order to ensure that the number of two-body collisions (that diverges approaching 
the critical temperature) remains small compared with the total number of test-particles 
(N co u « Ntp/2) at each time step. The factor F in eq.(4) gives the stimulated collisions 
induced by the statistics. For fermions the minus sign should be used instead of the plus 
sign; F = 1 leads to the classical Boltzmann collision integral. We note that similar Monte 



Carlo methods have been used for fermion |23].|24[ and boson [[25|] systems. 

The number of two-body collisions per unit time (collision rate) as a function of the 
temperature is given by: 
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n(p)n{p 2 ) [1 + n(p 3 )) [1 + n(p 4 )] (5) 
Expanding the B-E distribution as: n = (exp[(H — fi)/T] — 1) _1 = Y%Li ex P(j(f' ~ H)/T), 
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where H = 2— + U(r) with U the external (or mean-field) potential, we can carry out 



analytically all the integral in eq.(5), getting: 

cFcft (27r^,) 6 7rm V ; ^ i 2 j 2 n T J 

yi+j + fc (j+j + fc)(/x-y(f)) 
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Since we are considering homogeneous systems, we put U(r) = const. In fig. (1) we 
compare the analytical results for the normalized collision rate dN ^ n eq.(6) (dotted line) 
with the results obtained numerically with BMD (solid line). At temperatures T > 6 T c , the 
Bose-Einstein collision rate becomes equal to the classical Maxwell-Boltzmann one (dashed 
line). At T = 2.7 T c the B-E collision rate shows a minimum, and for T < 2.7 T c it 
increases rapidly deviating significantly from the classical Maxwell-Boltzmann trend. The 
fast increase is due to the stimulated collisions induced by the statistical factors (1 + n) in 
eq.(5). At high temperatures (classical regime), n « 1 and these factors are negligible, 
while, approaching the critical temperature, n » 1 and quantum statistic effects become 
dominant. To better illustrate this point, we show in fig. (2) the different contributions to 
the collision rate obtained from eq.(5). The long-dashed line corresponds to the term coming 
from n(p)n(p 2 ). This contribution decreases with temperature as it does the classical one 
(Maxwell-Boltzmann) in fig.(l). The dotted line corresponds to the term n(p)n(f>2)[n(p 3 ) + 
n(p 4 )], that shows the deviation from the classical trend. The biggest contribution to the 
collision rate comes from the term n{p)n{p2)n{pz)n{pi) (dashed line). 

The chemical potential ji is constrained by the normalization condition: 
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^expVumjl-W = (^)hoTl = 2.6(^)1 (7) 

3=1 



m T c 



where the critical temperature, obtained setting \x = 0, is T c = (2irh po/2.6 m)z. 

Eq.(7) shows that f = G[g where G is a universal function independent of the density 
of the system. This fact allows to calculate from fig.(l) the collision rate for an homogeneous 
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system with arbitrary density. In particular, for a given |r, we have: 



dN co u 



oc 



aTj' 2 (9) 



dfdt 

From the experimental point of view it is of particular importance to know the relax- 
ation time T re i ax , that is, the characteristic time for an anisotropic distribution to reach the 
equilibrium. On general grounds, it could be expected that r re i ax were proportional to the 
mean free collision time t co u. 



We define the relaxation time for a quadrupolar deformation as 26 



i = jxjjcoj m dp 

r re iax I n(p) P 2 {8) dp 
To evaluate eq.(10) numerically, we deform the initial equilibrium distribution as: 

p^p(l + aP 2 (6)) (11) 

where P 2 {9) = ^[3cos 2 (9) — 1] is the Legendre function. We note that for an arbitrary 
monopolar deformation (that implies only a scaling change in the modulus of the momenta 
p — > p(l + a)) the Boltzmann-Nordheim collision integral in eq.(l) becomes exactly zero 
(this holds also for fermions and classical particles). In this case the dynamics is entirely 
governed by the external (mean-field) potential. 

The quadrupolar relaxation time as a function of temperature is shown in Fig. (3) (solid 



line) for a = 0.5. We plot in the same figure the mean free collision time r< 



coll 



dN coH \-l 
dt > ' 



corresponding to the average time between two collisions for the system in the thermal 
equilibrium. These two characteristic times are proportional for T/T c > 3, but their behavior 
differs significantly at lower temperatures. For T/T c < 2.7 the collision time drops due to 
the increase in the collision rate (see fig.l) but, surprisingly enough, the relaxation time 
increases as: 

r rel <x(T-T c y 1/2 (2) 

Note that for a classical Maxwell-Boltzmann distribution the relaxation time goes as r^^ x oc 
y-i/ 2 . The reason for this behavior can be understood looking at the expression of the 



collision integral, eq.(l). The two terms involving the product of four distribution functions, 
F 4 = n(r, p)n(f, P2)n(f, p 3 )n(f, p^) cancel each other exactly for any distribution n(f,p,t). 
On the other hand, it is this term that, in eq.(5), gives the largest contribution to the collision 
rate at low temperatures. Moreover, the terms involving the product of three distributions, 
F 3 = n(r,fyn(r,p 2 )[n(r,p 3 ) +n(f, p 4 )], almost cancel each other for the deformation that we 
have considered. In fig. (4, a) we show separatedly the contributions of each term of eq.(l) to 
the relaxation of the initial quadrupolar deformation. The dot-dashed line indicates the time 
evolution of the quadrupolar momentum considering only the term F4 in the Monte Carlo 
simulation: as expected there is no contribution to the relaxation of the system. However, 
in fig.(4,b), we show how this term gives the largest contribution to the collision rate. The 
term F 3 gives also a small contribution in fig. (4, a) (dotted line), while it still gives a large 
contribution to the collision rate in fig.(4,b). Hence, we end up with the conclusion that the 
stimulated collisions do not contribute to the relaxation of the system. Only the spontaneous 
collision term (F — 1 in eq.(4)) affects the relaxation time, although this gives the smallest 
contribution to the collision rate (long-dashed line in fig.s(4,a,b)). 

As a consequence of these results, we remark that the simple criteria given in the litera- 
ture that state the transition between zero and first sound occurs when u t coU « 1, with ui 
the frequency of a collective mode, do not apply for a boson gas at T c < T < 2.7 T c , where 
T coii 0, but r re [ — > 00. Close to the critical temperature, the system should exhibit colli- 
sionless modes, with a transition to first sound at higher temperatures, when ur re i ax « 1. 

The predicted trend of r re i ax vs. temperature in fig. (3) can be tested experimentally. 
In particular, the damping of collective quadrupole motions, induced in the trapped gas, 
should rapidly decrease approaching T c according eq.(12). In effect the B-E distribution 
observed experimentally is non-homogeneous, being the boson gas trapped in an harmonic 
field. However we do expect that the conclusion of this work and the trend predicted by 
eq.(12) can still give a good approximation provided that the size of the system is much 
larger than the characteristic scale of inhomogeneity. 

We conclude remarking that the rapid increase of T re i ax for T — > T c questions the as- 
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sumption that the condensation observed experimentally takes place in equilibrated sys- 
tems. Understanding how the BEC occurs in systems far from equilibrium is one of the 
most interesting and challenging problems raised by the physics of trapped atoms. 

Part of this work has been done at the the Physics Department of the University of Illinois 
at Urbana-Champaign. T.L-A. thanks D.K.Campbell and A.S. thanks V.R.Pandharipande 
for the kind hospitality and support. 

Discussions with G.Baym, S. Fantoni, G.M.Kavoulakis, S.R.Shenoy and J.P.Wolfe are 
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FIGURES 

FIG. 1. Normalized collision rate as a function of the temperature, calculated numerically 
(solid line) and analitycally (dotted line). The dashed line shows the collision rate for a classical 
Maxwell-Boltzmann distribution. The particle density is p = 10 13 cto~ 3 and the collision cross 
section is a = 4.4 10~ 13 cm 2 . 

FIG. 2. Contributions to the collision rate from each term of the collision integral in BNE 
eq.(l), as explained in the text. 

FIG. 3. Comparison between the relaxation time (solid line) and the mean free collision time 
(dashed line) as a function of temperature. 

FIG. 4. (a) Quadrupole moment as a function of time for T = 100 nK. The different lines 
correspond to contributions coming from each term in the collision integral of BEN eq.(l), as 
explained in the text (see fig.2). (b) Corresponding collision rates. 
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